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Abstract 

Present evidences for neutrino masses and lepton flavour mixings allow to 
predict, in the Standard Model with light neutrinos, branching rates for the 
decays Z — > e/j,,fj,r, er of less than 10 -54 , while present experimental exclusion 
limits from LEP 1 are of order 10 -5 . The GigaZ option of the TESLA Linear 
Collider project will extend the sensitivity down to about 10~ 8 . We study in a 
systematic way some minimal extensions of the Standard Model and show that 
GigaZ might well be sensitive to the rates predicted from these scenarios. 
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1 Introduction 



Lepton flavour violation searches are as old as our knowledge about the existence of 
at least two different kinds of leptons: electron and muon. A prominent example of a 
lepton flavour violating (LFV) process is: 

H-^&y. (1.1) 

This reaction has not been observed so far, and the best experimental upper limit of 
its branching fraction is [|TJ: 

BR(/x -> e 7 ) = J}^ ^ C1 \ < 1.2 x KT 11 . (1.2) 
r(/i -> ^ez/ e ) 

At the Z factory LEP, searches for quite similar LFV processes, but this time 
directed to the Z boson, became possible: 

Z — > e/i, /xr, er. (1.3) 

The corresponding branching ratios are: 

BR( Z PA) = r '^f + ^» , (1.4) 

and the best direct limits (95% c.l.) are @: 

BR(Z -> e^) < 1.7 x 1(T 6 [3], (1.5) 
BR(Z — > e =F r ± ) < 9.8 x 1(T 6 1J, (1.6) 
BR(Z^ii T r ± ) < 1.2 x 1(T 5 (1.7) 

These (and many other) observational facts may be described with the concept of 
lepton flavour conservation (LFC) in neutral current reactions. In the Standard Model 
of electroweak interactions (SM) lepton flavour is exactly conserved. However, 

the model may be extended in such a way that virtual, LFC breaking corrections 
can appear. One mechanism relies on the assumption of neutrinos with finite masses 
and lepton mixing (from a non-diagonal mass matrix of the gauge symmetry eigen- 
states) IP-PII, leading to tiny rates for all the above processes caused by LFV one-loop 
effects. Historically, the i^SM — the Standard Model, enlarged with massive, mixing 
neutrinos — was the first theory allowing such predictions thanks to its renormalizabil- 
ity |12-14|. For the reaction fll.lj) and similar low-energy reactions like /i — > e~e + e~ 



or v\ — > V2I the first studies were reported in [|15|-|T7jj, and for the LFV Z decays flO] 
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The most general matrix element for the interaction of an on-shell vector boson 
with a fermionic current, as shown in Figure [I], may be described by four dimensionless 
form factors.0 At one-loop order, it is convenient to parameterize 

M = -^e'u f2 (p 2 )r p u h (- Pl ), (1.8) 



1 Soon later, related calculations were performed in the context of flavour non-diagonal quark 
production with a heavy virtual top quark exchange p0|-p2[ . 

2 For an off-shell vector boson two more form factors contribute. 



2 



fx (pi,mi) 

(Q)vwvQ 

^ /a (P2,m 2 ) 
Figure 1: T/ie effective LFV vertex. 

with a^y = g 2 /(47r), e being the boson polarization vector and 

T P = j p (fv - fAlb) + YF~ (ifM + fEl^Opv (1.9) 

Above, /v and stand for vector and axial-vector couplings and fu and f& for 
magnetic and electric dipole moments/transitions of equal/unlike final fermions. The 
form factors depend on the momentum transfer squared Q 2 = {P2 — Pi) 2 - For an on- 
shell photon, current conservation implies two additional conditions: (m.2 — rn\)fy = 
and (m 2 + mij/^ = 0. This means that LFV \i decays are exclusively due to dipole 
transitions, while for LFV Z decays all /y, /a, /at, /e are, in principle, non-zero. 

The general expressions for the branching ratios are:[] 

BR(/*-e7) = ^# (l/^ + i/il 2 ), (1.10) 



7r mf, 



BR(Z^) - ^ 



i/^r+i/.fp+^Vd/i/p+i/ii 2 ) 



Notice that while the muon total width is T M = aj^/(3847r) m^/M^, the Z width is 
T z ~ aiy/ c vF Mz- That is why BR(Z — ► £f^ 2 ) is naturally by an order of smaller 
than BR(/i — > ej). Furthermore, the M^/m 2 ^ enhancement of ( |1.1U| ) is compensated 
due to the chirality-flipping character of the dipole form factors, proportional to the 
fermion mass m^. 

The form factors are mo del- dependent. In the approximation of massless electrons 
(for fM — > ej) or massless external leptons (for Z — > l\l-f), there is only one independent 
form factor in each case. In the simplest assumption of n Dirac virtual neutrinos vi 
with masses m,j, the mixings factor out and one can write 



H -> e7 : 



-1*2 



fM = fl= ^^E V ^ V ^(A,; A Q ), (1.12) 



i=l 



1 n 

/m = /I = 0, (1.14) 



where V is the lepton-flavour mixing matrix and Vy/z are vertex functions, fully de- 



scribing the amplitudes. We have introduced the neutrino mass ratios A, = mf/My 



For the quark flavour-changing Z — > qiq2, multiply by a colour factor N c = 3. 
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and the virtuality of the Z boson Xq = Q 2 /Mw, that becomes \z = M\jM w on its 
mass shell.Q Owing to the unitarity of the mixing matrix, the amplitudes vanish for 
massless or degenerate virtual neutrinos, in exact analogy with the GIM cancellation 



in the quark sector 23 



We have strong evidence for neutrino masses of the order of some fractions of eV 
and large mixings |24],p5||. For small neutrino masses, a power-series expansion of the 
muon decay amplitude yields |~5l 



V y {Xi < 1; 0) « const 7 + A, + 0(A 2 ), (1.15) 

and similarly for the Z decayQ but with complex coefficients p!8|, p7| , |28[| : 

V z (\i < 1; X z ) ~ const z + (2.562 - 2.295 i)A< + 0(A 2 ). (1.16) 

The constant terms drop out after summing over the n generations of mixing neutrinos, 
but there survive contributions to the branching fractions proportional to the fourth 
power of the mass ratio mi/M w , for non-degenerate neutrinos, and thus unfortunately 
very small. Therefore, an observation of such LFV decays would be indicative to the 
existence of New Physics with a new, large mass scale involved. 

Consider now the hypothetical case of large neutrino masses. Neutrinos with large 



masses are accommodated by many extensions of the SM like grand unified theories [29 
or superstring- inspired models with an Eq symmetry [P0| . Heavy neutrinos are also 
well motivated by the seesaw mechanism PIH33I • From the exact expression of the 
LFV decays |1: 

V^AjjO) = Const 7 + $(A l ); (1.17) 
^ , . x(l — Qx + 3x 2 + 2x 3 — 6x 2 lnx) ,„ 

•w = u^y* • (L18) 



one obtains $(x ^> 1) — ► 2. In contrast, for the LFV Z decays |L9 



V z (Ai>l;A ) « Const z + ^ + 0(lnA 4 ). (1.19) 

Let us summarize the phenomenologically relevant differences between the LFV fi 
and Z decays: (i) the very different origin of the form factors intervening (dipoles in 
the /i case and mostly vector and axial- vector in the Z case); (ii) the 'typical size' of 
the rates due to the different powers of the coupling constant aw appearing in the 
branching fractions; and (iii) for fixed mixings, the Z branching ratio rises with virtual 
neutrino masses while the \x branching ratio reaches a plateau. 

In the rest of this work, we will concentrate on one-loop induced LFV Z decays. 
For these and other rare Z decays, the branching fractions are typically 

BR(Z^rare) ~ (— ) 2 ~C(10- 6 ). (1.20) 

4 The values M w = 80.41 GeV, M z = 91.187 GeV, c w = M w /M z , g = e/s w , and T z = 2.49 GeV 
will be taken throughout this work. 

5 This is in clear distinction to Eqn. (6) of | p6[ (with a logarithmic mass dependence), where from 
the recent neutrino data a prediction was derived to be BR(Z — > fJ^T ) ~ O(10~ 8 -j- 1CP 5 ). 
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There are many studies on such processes, in relation to e.g. CP violation 
neutral singlets W, 



, heavy 



j8|, supersymmetry p9| , f40| and superstrings |41| , |42fl or induced by a 
mixing with a heavy Z' |43j . See also the summary report of the LEP 1 Workshop |44| 
and the later study on the high luminosity LEP 1 project in particular [[4(J. The 
discovery reach of LEP 1 was indeed not very large, after comparing the experimental 
limits ( |1.5| )-( |ITT| ) with the order of magnitude of the potential effects ( |1.20| ). 

In a few years from now, a new high energy e + e~ Linear Collider could be con- 
structed. Interesting enough, with the GigaZ option of the TESLA Linear Collider 
project one may expect the production of about 10 9 Z bosons at resonance |j4?| . This 
huge rate, about a factor 1000 higher than the one at LEP 1, will make possible checks 



of the SM and its minimal supersymmetric extension MSSM at the two-loop level [48 



as well as searches for any kind of rare Z decays with unprecedented precision. A care- 
ful analysis [[4j| shows that in particular the LEP 1 discovery limits could be reduced 
to 



BR(Z 
BR(Z 



eV) 



< 2 x 10~ 9 , 

< k x 6.5 x 10" 

< k x 2.2 x 10" 



;i.21) 
1.22) 
1.23) 



with k = 0.2 1.0. This means one might have a chance of observation if the lepton 
mixings are not tiny and the masses of the neutrinos are at least of the order of the 
weak scale. Furthermore, in view of the expected sensitivities it might well be that the 
predictions are such that not only the asymptotic limit for large internal masses but 
an exact calculation of the effective vertex is needed: at least, it will be important to 
know where the large-mass limit fails. 

We perform a complete recalculation of the branching ratio (|1.4|) in presence of 
heavy Dirac or Majorana neutrinos and study the prospects for GigaZ in view of 
present, related experimental facts. We also compare to earlier studies and revise some 



of them. Many technical details of more pedagogical character may be found in |27j • In 
Section 2 the case of Dirac neutrinos is considered; Majorana neutrinos are treated in 
Section 3 and our conclusions are drawn in Section 4. The Appendix collects notations, 
conventions and useful expressions for the tensor integrals and the vertex functions as 
well as their low and large neutrino-mass limits. 



2 The LFV Z decays in the z/SM 

The simplest extension of the SM accounting for non-vanishing LFV Z decay rates 
consists of extending the particle content of the SM with three right-handed v singlets, 
thus forming three massive, mixing neutrino states a la Kobayashi-Maskawa. This 
is in conformity with compatible results from present solar, atmospheric, reactor and 
accelerator neutrino experiments. 

On basically the same footing one may also study the case of an additional se- 
quential, but heavy neutrino state. This case implies the existence of a heavy charged 
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lepton as well, in order to keep total lepton number L conserved.^ It is not a very 
favoured scenario but we consider it as a simple application. 

The final state charged leptons may be assumed massless. The amplitude is then 
purely left-handed and it is described by a single form factor, 

M = -J^ ; V(Q 2 )4u, 2 (p 2 )7p(l-75)u, 1 (-p 1 ). (2.1) 



Using the same vertex function Vz introduced in ( |1.13| ) one has: 



V Dir (Q 2 ) = ^V^V^A^Aq). (2.2) 



i=i 



In the 't Hooft Feynman gauge, the amplitude receives contributions from the set of 
diagrams of Figure ||: 

Vz{Xi]XQ) = v w (i) + v ww (i) + v^(i) + v^j,(i) + v w<t> (i) + v^(i). (2.3) 

The vertex diagrams Dl to D5 yield respectively: 

v w {i) = -(vi + ai) [X Q (Co + Cu + C la + C 23 ) - 2C 24 + 1] 

-(vi-ai) Xi Co, (2.4) 
v ww (i) = 2c 2 v (2/^)[Aq(C 11 + C 12 + C 23 )-6C 2 4 + 1], (2.5) 



A A 

v<t>(i) = -{vi + ai) y C - (vi - a,i) y 



Xq C 23 — 2C 24 + - 



(2.6) 



v H {i) = -(1 - 2s 2 w ) {2If) A, C 24 , (2.7) 
v w<p {i) = -2s 2 w (2P 3 L ) A, C . (2.8) 

The self-energy corrections to the external fermion lines DE contribute with: 

= ^(^ + a i -44a i )[(2 + A i )Bi + l]. (2.9) 

The definitions of weak neutral vector and axial-vector couplings are as usual: 

Vi = ^-23^ = ^(1-44^1), (2.10) 
a, = P 3 L , (2.11) 

and the dimensionless one- loop tensor integrals C , C , C a b, C a b and M 1 are given in 
Appendix ^, taking arguments Xi = Xj for the C functions. 

The form factor V describing the amplitude ( |2.1| ) is finite and no renormalization is 
needed, as expected because there is no tree-level coupling of a Z boson to two fermions 
of different flavour. Nonetheless, a non-trivial cancellation of infinities takes place, 
since C 24 , C 24 and Bi are UV-divergent. Actually, the vertex function Vz(Aj; Xq) is 
still infinite but has divergences independent of Xi, that makes possible the cancellation 
of the divergent terms in the amplitude, thanks to unitarity of the mixing matrix. 

6 A fourth generation of quarks is also needed to keep the theory anomaly free. 
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Dl: 




D4: vAAAAX,' 



DS: 



Figure 2: Feynman diagrams for the lepton- flavour changing Z decay. In the case 
of virtual, ordinary Dirac neutrinos, the ZuiUj vertices in Dl and D3 are flavour- 
diagonal. The analogous quark-flavour changing process can be obtained by replacing 
Ik by down-quarks and Vi by up-quarks. 

2.1 Contributions from light neutrinos 



Disregarding the controversial results of the LSND accelerator experiment, all neutrino 
experiments are compatible with the oscillation of three neutrino species. We will 
now estimate the LFV branching ratios under the assumption that there are three 
generations of light neutrino flavours and that their mixing is given by the unitary 
mixing matrix V constrained by the experiments. The mixing is described by three 
angles $12, $13, $23, and one CP-violating phase 5 as in the quark CKM case.Q 

A global analysis of atmospheric neutrino data favours i/„ — v T oscillations 



Am 
sin 2 2-d 



atm 



' 28 



(1 6) x l(T 3 eV 2 



at i n 



sin 2 2^23 ^ 0.8 4- 1.0. 



(2.12) 
(2.13) 



The solar neutrino deficit is compatible with v e — oscillations f50f h 



sin 2 2dp, 



Am, - 



i\ 2 ~ 1(T 10 4 10" 
sin 2 2$i2 = free. 



(2.14) 
(2.15) 



There are solutions for vacuum and matter oscillations compatible with a wide range 
of masses and mixing angles, although the large mixing angle solution LMA with 



7 Oscillation experiments cannot distinguish between the Dirac or Majorana character of the neu- 
trinos. If they happen to be Majorana particles, two additional CP-violating 'Majorana' phases a, j3 
are needed since for strictly neutral particles less phase factors may be 'eaten' by redefining complex 
fermion fields. They are set here to a = j3 = 0. 
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maximal mass splitting seems favoured. From reactor searches, there are no hints of 
v e — v T oscillations which allows us to assume 



sin 2 2d 



13 



0. 



(2.16) 



Taking this information into the standard parameterization for the mixing matrix 
one has 



V 



Cl2 
1 



Sl2 





1 



y^ 12 7f 12 si 

73 Sl2 ~72 Cl2 75' 



(2.17) 



Using the unitarity of V and t\ ^ £ 2 , 

3 



BR(Z 



a 3 w M z 
\Wc 2 w Y z 



E v ^ v u^( A ^)-^(Mz) 



i=i 



[2.18) 



Performing a well justified low neutrino mass expansion of the vertex function (see 
Appendix |0|), one finds 0^7j|: 



V z {X i ,X z )-V z {0,Xz) 



ai Xi + 0{X\ 
2.5623 - 2.2950 i. 



(2.19) 
(2.20) 



1 ^) goes as m\ for low neutrino masses. This behaviour is shown 



Therefore BR(Z - 

in Figure It is valid over a large mass range until about m 8 « 30 GeV, i.e. just 
below the Z mass. 

Taking now into account the phenomenological squared mass differences Xij = 
Amfj/M^r and the mixing angles ( 2.12 )— ( ^.16 ), one can determine the finite expec- 
tation: 



BR(Z -> /x T r ± ) ~3x 10" 6 x |sf 2 Ai 2 - A 23 | 
and the upper limit: 

BR(Z -> e T / u ± ) ~ BR(Z 



;i -=-30) x 10 



-55 



6 x 10~ 6 x c\ 2 s\ 2 \\ 2 < 4 x 10" 60 . 



(2.21) 



(2.22) 



These extremely small rates are far beyond experimental verification. This justifies 
taking the light neutrino sector as massless in the following sections where we discuss 
extensions providing larger rates. 



2.2 Contributions from one heavy ordinary Dirac neutrino 

Assume the neutrino of generation N to be the only heavy one, mixing with a light 
sector with negligible masses. Then, using again the unitarity of the mixing matrix: 

BR(Z^ifit) = 1Q »Vf r \^ Y tN\ 2 \Vz(X N ;X z )-V z (0;X z )\ 2 . (2.23) 

iy/7r c w l z 
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BR(Z -> if if) with unit mix. 



v u Q 2 = M 2 



0.01 




10" 



10" 



10" 



10" 



10" 



10" 



10" 



10" 



BR(Z -> with unit mix - / 



i = i/ is Q 2 = M 2 





i = Ui, Q 2 = M\ ft 
i = Vijm, Q 2 = // ,. 



100 
mi [GeV] 



1000 



Figure 3: Contribution of one neutrino generation i to the LFV Z 



1 *-2 



decays 



for ordinary Dirac neutrinos in the small and large neutrino mass regions, and the 
analogous quark case. The mixing factor has been set to V^jV* ■ = 1. 



In the large Dirac neutrino mass limit, the following approximation works well (see 
Appendix |A.2| ): 



Vz(\; Xz) 



with 



-4ci / A e + A + I 3 ^-(1 - 2s 2 



w , 



In A 



(-66 - X z + 9Qs 2 w + 5s 2 w X z ) 



y/l/\ z - 1/4. 



(2.24) 



+ i (-8 + 2A Z - 324, - 4s 2 w X z ) y arctan (J- J + O (^) 



(2.25) 



The vertex function contains a constant term proportional to A e = 2/(4 — c?) + 7+ln47r, 
divergent in d = 4 dimensions. This term drops out in the physical amplitude, as ex- 
pected, since the unitarity of the mixing matrix demands the subtraction of V z (0; Xq), 
with identical divergences. Its expression can be found in Appendix A.l| . For an 
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on-shell Z, 



Vz(\n, \z) - V z (0- X z ) = -[X N + 2.88 In A,v - (6.99 + 2.11 i)} + 0(\n X N /X N ). (2.26) 



The exact results are depicted in Figure |3|, where the simpler calculation with Q 2 = 
is also displayed. We find agreement with earlier calculations [T§[], also for quark 



flavour- changing Z decays pl|, 22 1. Of course, the results for Q 2 = are a good 
approximation only when m 2 N 3> M§. 

For a study of the size of the branching ratios, the knowledge of the light-heavy 
mixing elements involved in ( |2.23p is crucial. Their values do not only influence poten- 



tial LFV processes but also flavour-diagonal ones. Using a general formalism developed 
in one can exploit measurements of flavour diagonal processes (checks of lepton 
universality and CKM unitarity, Z boson invisible width, etc.) ]54|J55[] to obtain indirect 



experimental bounds on such light-heavy mixings pql , defined as 



[2.27) 



The most recent indirect bounds [57| : 

(2.28) 
(2.29) 
(2.30) 

are only improved by the impressive accuracy of the direct searches for LFV processes 
involving the first two lepton generations. In fact, for heavy enough neutrinos one can 
rewrite (|1.10|) : 



< 


< 


0.012, 


< 


< 


0.0096 


< 


< 


0.016, 



3a 



BR(/x - ery) « ^-s 2 v / 



and from ( |1.2| ) a stringent 'mass-independent' limit can be extracted |58 

s 2 s 2 < 1.4 x 10~ 8 . 



(2.31) 



(2.32) 



The limit above sends the Z — > e/x process beyond any experimental reach, even if 
the neutrinos are very heavy. At this point, it is important to realize that, although 
the branching fractions for large neutrino masses grow as 

BR(Z^44)ocs^^ 2 m^, (2.33) 

there is a 'natural' upper value for the neutrino mass determined by the perturbative 
unitarity condition on the heavy Dirac neutrino decay width |59|-|63| . 



10 



that implies 



m 



N 



< 



4M^ 
aw 



£ 

.fc=i 



(4.4 TeV) 2 . 



(2.35) 



In other words, expression ( |2.34|) shows that the unacceptable large-mass behaviour of 
the amplitudes (oc m 2 N ) is actually cured when a sensible light-heavy mixing (at most 
oc mjj 2 ) is taken into account |64].|65|. 



l(T 3 
10- 4 
l(T 5 
1(T 6 
10- 7 
10- 8 
10" 9 
10- 10 
10- 11 

lO" 12 
lO" 13 



BR(Z -5- pFr* 
Dirac ordinary 



Dirac singlet 

m N .Jm Nl = 2 
m N .Jm Nl = 3 
mjv 2 / m iVi = 10 




100 

mjvj [GeV] 



10000 



Figure 4: Upper limit of BR(Z — > /U =F r ± ) assuming a light neutrino sector mixing with: 

(i) one heavy ordinary (thick- solid) or singlet (thin-solid) Dirac neutrino of mass m^; 

(ii) two heavy right-handed singlet Majorana neutrinos (dashed lines) with masses 
and mN 2 ■ 

For illustration, we show the less constrained results for Z — > /i =F r ± in Figure § 
(thick-solid line), assuming the present (indirect) upper bounds on the corresponding 
mixings ( p.29|) and ( |2.30| ): given the heavy neutrino mass(es), the branching fractions 
cannot exceed either the curves or the collider exclusion limit . Regions below the curves 
correspond to mixings smaller than the upper bounds ( |2.28|) to (|2.30| ). Masses beyond 



the end points are acceptable only if the mixings are smaller than the upper bounds. 
Of course, since there are at least two unknowns, a neutrino mass and a combination 
of mixings, the LFV Z decays cannot improve the bounds on the mixings without 
assuming a value for the heavy neutrino mass(es). This is in contrast to the LFV [i 
decays for sufficiently heavy neutrinos fl2.31| ). 
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3 The LFV Z decays in the z/SM with right-handed 
Majorana singlets 



Let us now consider the case when the heavy neutrinos are Majorana particles. Actu- 
ally this a very interesting possibility since such states belong to the particle content 
of most GUT theories, like SO(10). Furthermore, they may participate in the seesaw 
mechanism, that explains the smallness of the observed neutrino masses by introducing 



a general Majorana neutrino |66| mass matrix, incorporating ordinary Dirac mass terms 
Trip, of a size typical to the charged lepton sector, and lepton-number violating Majo- 
rana mass terms at a higher scale M R 3> mjj. Majorana mass terms MrV r Vr + h.c, 
with vr being right-handed singlets under the SM group, are gauge invariant, but vio- 
late lepton number by two units. The physical states after diagonalization of the mass 
matrix are, respectively, light and heavy Majorana neutrinos with masses 

m v « m 2 D /M R , m N « M R > m v . (3.1) 

If there is only one generation of heavy neutrinos, the light-heavy mixings are fixed 
to be very small, 



m D /M R « ^Jm v /m N , (3.2) 



leading to unobservable LFV effects. 

But this is not the case when one includes several right-handed Majorana neutrinos 



with inter-generation mixings ||56|,p7|,|68|. We will focus on the most conservative case 
of two heavy right-handed singlets. 



3.1 LFV with Majorana neutrinos 

Let us consider n G generations of charged leptons (Dirac fermions), whose left-handed 
components (£° = &Li ■ • •) belong to the same isodoublet as n G left-handed 

neutrinos (v® = u e , u T , . . . ) and, in addition, n R right-handed neutrino singlets. 
The interaction eigenstates are a mixture of physical states given by p6| , |69| ,[70| 

t° Li = E U &^> ( 3 - 3 ) 

riG+ n R 

"li = E V M> ( 3 - 4 ) 
3=1 

where v=7] v c are n G + n R Majorana fields (i.e. self-conjugate up to a phase rj). 

In the charged- current interactions, one must replace the leptonic mixing matrix V 
by its generalized version, the rectangular n G X (tiq x n R ) matrix B, 

n G 
k=l 
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Therefore, in the physical basis, 



Ccc = ^WAH^l + h.c. 



(3.6) 



where P R%L = |(1 ± 75). 

But the main feature distinguishing Dirac and Majorana cases is the existence of 
non-diagonal ZviVj vertices (flavour-changing neutral current), coupling both left- and 
right-handed components of the Majorana mass eigenstates to the Z boson, 



-C 



9 „ or 



NC 



2c w 



Z.K^PlvI-vI'^Pr^ 



9 ry — 



2cw 



Z^C^Pl - CWPr 



v 



.1 * 



(3.7) 



where = Gv\ is the charge-conjugate of z/°, which is right-handed, and 



n G 



C ij = U fci Ufc i' = n G + n R ), 



(3.1 



k=l 



a quadratic (uq + n R ) 2 matrix. Such flavour-changing NC vertices appear in graphs 
Dl and D3 of Figure |2] where Majorana neutrinos couple directly to the Z, and a W 
or a Goldstone boson (ft is exchanged: 



-C ij [X Q (C + C 11 + C 12 + C 23 

2 



2C 24 + 1] + C* ij y/X~X~ Co, (3.9) 



— 2C24 + - 



(3.10) 



The other diagrams remain unchanged compared to the Dirac case and the resulting 
form factor reads:^ 

V Maj (g 2 ) = J2 Bii&tvVzihj), (3.11) 
».i=i 

V z (i,j) = v w (ij) +v^(i,j) +v ww (i) +v H (i) +v W(f> (i) +vj:(i). (3.12) 



We have used the Feynman rules in W% - V% 
Majorana particles. 

It turns out convenient to cast (|3.12j) as 



to properly handle interactions involving 



Vz(Ai, Xj) = 6 t3 F(X l ) + C t3 G(X t ,X J ) + C* J ^X~X j H(X l ,X J ). (3.13) 



8 We have compared our formulae with Eqn. (Bl) of fl7lfl and found disagreement, in particular 
the appearance of a tensor integral C22 at several instances. C22 is UV-finite and has no numerical 
impact on the amplitudes for large neutrino masses, where we find full agreement. However, the 
rearrangement of (3.14) leads to the well established Dirac vertices only when using our expressions. 
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The Dirac vertex function ( |2.3j ) is then 

V z (Xi) = F(X i )+G{X ii X i ). 



(3.14) 



The form factor ( |3.11| ) is UV-finite, but the vertex function Vz{Xi,Xj) is not. The 
divergences are such that they exactly cancel due to unitarity relations among the 
mixing matrix elements of B and C [ 27 , p6(1 . The same relations allow to write VMaj in 
terms of only the heavy sector, assuming the light sector being massless: 



V Maj (Q 2 



n R 

B^B^. 

j,i=l 

x {5 NiNj [ F(X Ni ) - F(0) + G(X Ni , 0) + G(0, X Ni ) - 2G(0, 0)] 

+ C^. [ G(X Ni ,X N .) - G(X Ni , 0) - G{0, X N .) + G(0, 0)] 

+ CjViiv,- y A^. H (X n . , X Nj ) [ . (3. 15) 



3.2 The z/SM with two heavy Majorana singlets 

In the simple case of tir = 2 heavy right-handed singlet neutrinos Ni and N2, mixing 
with a massless sector, the B and C matrices are fully determined by the ratio of the 
two physical heavy masses squared r = m 2 N<2 / m 2 Nl and the light-heavy mixings s 2 , , here 



(3.16) 



Their explicit values are |56 



C 



W1JV2 



B^ 

B^AT 2 
G>N 2 N 2 

Cn 2 n 1 



,1/4 



i 

Vl + rV2 

r l/2 



1 + r 1 / 
1 



1 + rV2 
1 + rV2 



1/2 

E 



E 



(3.17) 
(3.18) 

(3.19) 
(3.20) 
(3.21) 



The mass ratio r is a free parameter and the light-heavy mixings are constrained by 
present experiments as shown in Section [2.2| . Upper values for the branching ratios of 
Z — ► / Li T r ± , obtained from the experimental bounds given the heavy masses m Nl , mjv 2 , 
are also displayed in Figure [|. 

The case of two equal-mass Majorana neutrinos = m^ 2 is equivalent to one 
heavy singlet Dirac neutrino,^ and it approaches rapidly the ordinary Dirac case for 

9 In fact, two equal-mass Majorana neutrinos with opposite CP parities form a Dirac neutrino. 
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small masses. This phenomenon is just another example of the "practical Dirac- 
Majorana confusion theorem" |74[ (see also the recent discussion in |75|,[7C| and ref- 
erences therein). If both masses and rriN 2 are small, the amplitude goes as 
V^AiVi = rn^m^/M^r with the same global factor a\ as in the ordinary Dirac case 
(|2.20| ). This can been seen in Figure (| not far below the Z peak, where the branching 
ratios grow with \ 2 N and scale with the ratio of the two neutrino masses squared. 

If one of the neutrinos has the mass of the Z boson, the imaginary parts of the 
amplitudes VDir,Maj (coming from the subtraction(s) at A at = 0) dominates, both for 
the Dirac and the Majorana cases. This happens since the real parts are slowly varying 
for M/v < Mz, while the imaginary parts vanish for Mz < + M^i. Further, since 
these imaginary parts necessarily come from accounting the subtractions of the zero 
mass limits, they are independent of the value of r. This results in common values of 
the branching ratios for rriN 1 = Mz for any value of rriN 2 . Nevertheless, the subtraction 
of the light sector implied by the unitarity constraints is not the same for the cases of 
a heavy ordinary Dirac neutrino and heavy Majorana singlets. One finds explicitly 

= -1.0524, ^_J^1 = -2.0653. (3.22) 



The expansion of the form factor (|3.15 ) in the large neutrino mass limit A^ ^> 1 
at fixed r, leads to (see Appendix |A.2| 



VMaj (Q z 




(3.23) 



that agrees with |56j for the unphysical value Aq = 0. The constant in front of the 
leading term coincides for r = 1 with the ordinary Dirac case, except for an extra 
damping factor J2i s l e > that makes the Dirac singlet case in particular, and the Majo- 
rana case in general, more sensitive to the present bounds on the light-heavy mixings. 
The constant in front of the In A term, subleading but not so much mixing-suppressed, 



is identical to the one in the ordinary Dirac case ( |2.24 ) 



We have cut again the curves at the perturbative unitarity mass limits. Due to the 
different number of degrees of freedom of the Majorana particles, the condition on the 



heavy neutrino width is this time [52 



2 x ^rn% |B Wl | 2 < \m Ni , (3.24) 
resulting in the mass limits 



m% = -m 2 N2 < \±ll!l x (3.1 TeV) 2 . (3.25) 



r 



We see from the figure that GigaZ has a discovery potential, preferentially in the 
large neutrino mass region, if the light heavy-mixings are not much below the present 
upper limits. Due to the different coupling structure, the simple sequential Dirac 
neutrino case does not constitute a limiting case for large masses. 
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4 Concluding remarks 



The sensitivity of the GigaZ mode of the future TESLA linear collider to rare, lepton- 
flavour violating Z decays has been studied. We have determined the full one-loop 
expectations for the direct lepton-flavour changing process Z — > £i£ 2 with virtual Dirac 
or Majorana neutrinos. This is an interesting theoretical issue in view of the evidences 
for tiny neutrino masses from astrophysics, which might be also indicative for the 
existence of heavy neutrinos in some Grand Unifying Theory. Both the exact analytical 
form and the large and small neutrino mass limits of the branching ratios are given, 
thereby cross-checking the existing literature. From our numerical studies, taking into 
account the present experimental results, we conclude that: (i) the contributions from 
the observed light neutrino sector are far from experimental verification (BR < 10 -54 ); 
(ii) the GigaZ mode of the future TESLA linear collider, sensitive down to about BR 
~ 1CT 8 , might well have a chance to produce such processes, if heavy neutrinos exist 
in Nature and if they mix with the light ones in a sizeable way. Finally, we have shown 
that we could gain from observation of the LFV Z decays alternative informations 
compared to the LFV \i decays. 
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A Tensor integrals and vertex functions 



We have introduced dimensionless two- and three-point one-loop functions: 

Bi(Ai) = B^O-mlMl), 
C..(Ai) = M^C(0,Q 2 ,0;m2,M^,M^), 
C.XX^Xj) = Mf v C.X0,Q 2 i^Mt v ,mlm 2 J ), 



(A.l) 
(A.2) 
(A.3) 



from the usual loop integrals 1 14], [77 1 with the tensor decomposition (Minkowski metric): 



pfBu (A.4) 
rfC n +j$C 12 , (A.5) 
fxP\C 2X + f 2 p\C 22 + (pft£ + pMCaa + <T A&CW 

(A.6) 



The tensor integrals are numerically evaluated with the computer program LoopTools 



[ 78] , based on FF ]75 |. All the numerical results for the Dirac case have been carefully 



checked against an older approach described in |L9] 
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The following definitions of the integrals in d dimensions are useful: 

li(Ai) = 



A e f 1 

— - + / dx xln[(l — \i)x + Xi — ie], 
2 io 



A e 1 f 1 f x 
X _ 2 ' dx ' dy lnZ? « w ' 



o jo 



f dx f dy [l,-y,y(l-x)) ——, 

JO JO UijW 



C 2 4(Aj, Xj) 

Co ,11,23 (Xi, Xj) 

with Xi = m 2 jM^, A t = 2/(4 - d) - 7 + ln4?r and 

Dijw = ^zxy + (1 - Xj)x + [-X z + (Xi - l)]y + Xj - ie. 
To get the barred tensor integrals C, replace D^w by: 

D W Wi = X z xy - (1 - Xi)x + [-X z - (Xi - l)]y + 1 - ie. 



(A.7) 
(A.8) 
(A.9) 

(A.10) 
(A.ll) 



The functions Bi, C 2 4 and C24 are UV-divergent but the physical amplitudes are 
finite. 



A.l Light neutrino mass expansions 

Let us first list the value of the necessary tensor integrals for massless neutrinos and 
Aq + @: 



Bi(0 
C (0,0 
Cu(0,0 

c 12 (o,o 
c 23 (o,o 

c 24 (o,o 

C (0 
C n (0 

C 12 (0 
C 23 (0 

c 24 (o 



_Ae_ 1 

2 4' 
-c , 

-T- (co - 1+lnAg -m) 
Cu(0,0), 

1 



(A Q + 2)c -^-2 + 2(lnA c 



m 



(A.12) 
(A.13) 
(A.14) 
(A.15) 
(A.16) 



"Ob 



T + 4A~ h2(AQ + 1)Co + 3Aq + 2 ~ (Aq + 2)(ln Xq ~ (A ' 17) 

(A.18) 

(A.19) 
(A.20) 
(A.21) 



Cu(0) 



(c -5 + 1), 



A, 



4 2Ar 



c - B + 1 - -A Q + 7rA Q y - 2A Q y arctan(2y) 



(A.22) 
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with 



Aq c 
Aq c 



6 U + A c 



ln 2 (l + A Q ) + 7rln(l + A Q ) z, 



(A.23) 



7T 



Li 2 (l - A Q ) + 2KeLi 2 



(Ac 



A, 



1 + XqV i 



A 



- 23?eLi 2 ( 1 - ^ - X Q y i ) , 



B 



2?7 



arctan(2y) + arctan ( — — - — 2y 

3 — Xn 



(A.24) 
(A.25) 



After expanding the tensor integrals for small neutrino masses (see Appendix D.2 
of ||27||), the vertex function for the case of a light Dirac neutrino reads: 

(A.26) 



V^(A<1;Aq^0) = V z (0; X Q ) + ai X + 0(\ 2 ) 
where the terms proportional to A In A have cancelled out and 



V z {0] A 



Q) 



O _|_ Q \ I 

^2c 2 w A e + -4^(mA -7T i) - W (7X 2 Q + 14A Q - 



2Ar 



;i + a q ) 2 



Co 



+ + 2A Q )(c -B) + ^-[iry- 2y arctan(2y)]. 



(A.27) 



Only the functions C . develop imaginary parts, and only for Aq > 4A«. At the Z peak 
the numerical result is: 



V z {0;X z ) = -2c 2 w A t + 1.2584+ 1.0524 i. 



(A.28) 



The linear term in the expansion (|A.26 ) has the coefficient [18, 27], 28 



2 1 2 

ai(X z ) = -— (1 + X z )c + 7TT2"( 4 A| - 5A Z - 6)c - — (\nX z - n i) 
\ z Z\ z \ z 

+ ^r(25A| - 38A Z - 24) + JL(2 - X z )ny 
+ T2"(A| + 7X Z + 6)y arctan(2y) 
+ t-o"(3Az + 2)y arctan ( —2?/ 



X 2 

2.5623- 2.2950 i. 



3- A; 



(A.29) 
(A.30) 



The behaviour of ( |A.26j ) is in contrast to the case Xq = for which |I9~, 27\, 52 

V Z (X < 1; X Q = 0) = i(-4c^A e + 6A + 2A In A) + 0(A 2 ). (A.31) 
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A. 2 Heavy neutrino mass expansions 

The limits of the necessary tensor integrals and the vertex function in the Dirac case 



for large neutrino masses can be found in Appendix D of ||27|| . We collect below the 
large mass expansions of the tensor integrals that are also needed for the Majorana 
case, namely one or two identical neutrinos running in the loop: 



Cu(Ai 



Bi(At) 
C (A,) 
C 12 (A,) 
C 23 (A,) 
C 24 (A 4 ) 



In \j 

Xi 

lln Xi 
1 In A, 



1 — 4y arctan I — 



6 A, 



+ 0(1/X 2 
+ 0(1/X 



I ) ■ 



1 + 0(1/*?), (A.33) 
(A.34) 
(A.35) 



In Aj 
12A" 



+ 



-30 + 5A Q + 24(4 - X Q )y arctan ( — 



72A, ; 



2 ) 

l ) ! 



Co(Ai, Xi) 
Cn(Aj, Xi) = Ci 2 (Ai, Xi) 
C 2 3(Aj, Aj) 
C 2 4(Aj, Aj) 

with 



+0(1/X 
1 In A; 



1 

4A~ 



12A? 



0(1/X 



i ) ■ 



18A,; 



o(i/x 2 ), 
+ o(i/x 



2 ) 



A e 1 , , 1 , 

f -4^A i + - + (-9 + A ( 



36A, 



y= ^I/Aq-1/4. 



(A.36) 
(A.37) 

(A.38) 
(A.39) 

+ 0(1A 2 ), (A.40) 



(A.41) 



Substituting the expressions above in ( |2.3| ) one gets the Dirac vertex function of fl2.24|) . 

Besides, we need some additional expansions for two Majorana fermions with dif- 
ferent large masses masses Aj ^ Xj, 



Co(Aj, Aj) 



Xi — X-i 



A, 



+ 



X Q 
(A, - Xj) 2 



In Xi 



2Ai \ , x 1 A 2A, 



A,' lnA ^' 



+ C(l/A 3 



(A.42) 
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and, to a lower accuracy in the expansion parameters: 
Cn(Aj, Xj) 



Cu(Xi, Xj) 
^23(^4, Xj) 

C2/t(Ai, Xj) 



11 + lnAi , 

1 — - 

2 Aj — Aj 

Cn(Aj, Aj), 

1 



A, ' luA ' -AilnAi +0(l/A 2 ), 



2(Ai - x 3 y 



6(Ai - a,; 



2A A 

(Aj + Xj) - — -4- (In A*- In A, 



Ai Aj 



+ 0(1/A 2 



(A.43) 
(A.44) 
(A.45) 



3 l(A < + l)lnA<-(Aj + l)lnA j 
8 _ 4 Ai - A~ 



^C a8 (A i ,A i ) + 0(l/A 2 ).(A.46) 



Actually, Cn and C12 are irrelevant for large neutrino masses. 

Finally, in ( 3.1 5[) we need loop integrals where one neutrino mass is large and the 
other one vanishes. They are all irrelevant except C24 in this limit, but we show their 
expansions for completeness: 



C (A,0) 
C 12 (0,A) = C n (A,0) 
C n (0,A) = C 12 (A,0) 
C 23 (A,0) 
C 24 (A,0) 



In A In A 



A 



A 2 



+ 



^5.4 nc\ /> :1 
2A 2 



2^ + (3 + Aq) 6A^"^ "" ,V;i 



+ 0(1/X i 

, ■ on. a- 



In A 

6A 2 " 



12A 2 

^ + -- — -^- + 0(1/A 2 ), 
2 8 4 12A w h 



+ C»(1/A 3 



(A.47) 
(A.48) 
(A.49) 
(A.50) 
(A.51) 



and 



Co ,23,24(0, A) — Cq,23, 24(A, 0). 



(A.52) 
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